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Abstract 

A new method for elimination of the ’spurious’ states comprising requirement 
of translational invariance and simplicity of the enumeration scheme of antisym¬ 
metric A-particle states has been developed. The method presented enables one 
to project off the space of purely spurious antisymmetric A-particle oscillator 
functions with singled out dependence on intrinsic coordinates of two last par¬ 
ticles onto the subspace of the nonspurious motion. The method is based on 
constructing the matrix elements of the centre-of-mass Hamiltonian and spectral 
decomposition of the corresponding real symmetric matrix. The obtained coef¬ 
ficients of the expansion of the antisymmetric A-particle oscillator functions in 
terms of the ones with singled out dependence on intrinsic coordinates of two 
last particles and with eliminated spurious states are required for calculation of 
two-particle oscillator intrinsic density matrices. Software tools based on a root 
rational fraction expression of numbers for the theoretical formulation presented 
have been developed. 

1 Introduction 

The translationally invariant shell model, originally proposed by Kretzschmar |^, 0, 
offers a very sophisticated method for characterization of antisymmetric states and 
producing the intrinsic coefficients of fractional parentage (ICFP’s). A description 
of the many-particle wave functions in terms of relative coordinates turns out to be 
unfeasible because of the antisymmetrization requirements of the Pauli principle. As 
a consequence the ICFP’s are available only in the E min approach and for the excited 
states with the maximum orbital momentum |^, 

More promising is to establish and use the relation of the translationally invariant 
shell model (for an arbitrary excitation) and the ordinary oscillator shell model. The 
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relation of this type makes it possible to obtain the results that obey the principle of 
translational invariancy by means of calculations within the framework of considerably 
simpler oscillator shell model. 

However, the ordinary oscillator shell model is inherent with the so-called ’spurious’ 
states corresponding to the centre-of-mass motion of the nucleus |^. The ordinary 
oscillator shell model Hamiltonian is not invariant with respect to translations. Hence, 
the wave functions under calculation may contain nonphysical components, which are 
called spurious. The latter describe excitations of the centre-of-mass motion of the 
nucleus that are introduced with the single-particle potentials in the laboratory system. 
These excitations should be eliminated, since only the intrinsic excitations corresponds 
to the genuine excitations of the physical system that are observed. 

The description of intrinsic properties of atomic nuclei may be simplihed using the 
two-particle intrinsic density matrices instead of wave functions |^. However, according 
the procedure proposed in Ref. 0 , for calculation of the two-particle intrinsic density 
matrix for an arbitrary number of oscillator quanta, the expansion coefficients for the 
product of the centre-of-mass ground state function and the intrinsic wave function in 
terms of antisymmetric H-particle functions are needed. 

In the present work we show that the coefficients of this type may be obtained by 
generating purely spurious antisymmetric H-particle oscillator functions with singled 
out dependence on intrinsic coordinates of two last particles and, then, by projecting 
them off onto the subspace of the nonspurious motion. This may be accomplished 
by diagonalizing the centre-of-mass Hamiltonian in the basis set of antisymmetric A- 
particle oscillator functions with singled out dependence on intrinsic coordinates of 
two last particles and choosing the subspace of its eigenvectors corresponding to the 
minimal eigenvalue equal to |. This procedure can be successful only when the non- 
truncated conhguration space is taken into account. It implies that a large number 
of antisymmetric H-particIe states has to be used. Such large scale calculations can 
be carried out by using the algorithms presented in Ref. |^, for enumeration and 
formation of antisymmetric H-particle states. The proposed method for elimination of 
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’spurious’ states is the first one which enables practically production of two-particle 
oscillator intrinsic density matrices for an arbitrary excitation. 

2 The centre-of-mass Hamiltonian matrix 

The functions from the subspace of the nonspurious motion are the eigenfunctions of 
the centre-of-mass Hamiltonian. Thus, projection onto the subspace of the nonspu¬ 
rious motion implies such superposition of the oscillator shell model functions which 
diagonalize the centre-of-mass Hamiltonian matrix. So, at hrst the matrix elements 
of the centre-of-mass Hamiltonian between the osillator shell model functions (in the 
conhguration space) should be obtained, at second, the eigenvalue problem should be 
solved, and at last the submatrix corresponding to the zero excitation energy should 
be extracted. 

The Hamiltonian of the many-particle oscillator shell model may be divided in two 
parts: the centre-of-mass Hamiltonian and the intrinsic Hamiltonian Tfintr 

^ ~ -^c.m. + -^intr. ■ (1) 

The centre-of-mass Hamiltonian describes the motion of the centre of mass 

of the H-particle system and is of the form 

= f {-SvUpfld ■ P) 

Here R = ^ Hf=i denotes the centre-of-mass position vector, Vr = J2t=i ^ 

related with the usual oscillator length via b = where m denotes the nucleon 

mass. In the one-particle coordinates the Hamiltonian (2) can be written as 

= is I (e V,) + ^ (e I . (3) 

By using b as our unit of length and rewriting Eq. (0) in terms of dimensionless 
quantities = Ti/b and Vp) = 6Vj, one obtains 


3 



( 4 ) 


-^c.m. 2A X! { ^(0^0) + ^(0^0)} • 

*j=i 


In order to separate mixed terms one can use the identity valid for an arbitrary entries 
ttij symmetrical with respect to the index permutation aij = aji 


A A 

op' = 4" ~ + 

ij=l i,j=l 

i<j 


■ 


( 5 ) 


This leads to representation of -ffc.m. terms of single-particle and two-particle opera¬ 
tors 


TT _ hu) 2 

-^c.m. — 2 A(A-l) 


E {^1 [-v?p + .1) 


^ +A| 


~^u) + yu). 


*j=i 

i<j 


[-(V,., - Vu,)'^ + (4 - 4,)]} 


( 6 ) 


Now, the matrix element of the centre-of-mass Hamiltonian between the oscillator shell 
model functions assumes the form 


{EKAjnT\- 


H. 


c.m. 

huj 


lEK'A'JHT) 


E 


y] {{EKAJETy, {{elj)A-i, {elj)A)rT"\\EKAjnT) 


(EKAjnT) 

j/fj'/f (^0) A _ 15 (e/ j) ^ 


X {(EKAjnry, ((eij)Pi, (elj)'AJ"T"\\EK'A'jnT) 

X {((el])A-u (dj)A)rT"\ {A\ [-V 4 „ + 4 _,|] + Ai [-V^ + 4 , 

-(/I - i)i [-(V,x-1) - v,x,)= + (ma-,) - »(X))T l(W)Vi. {eijurr") . (i) 

Here the usual shell model fractional parentage expansion is used to separate out two 
last particles in the list, and the notation follows that introduced in Ref. |^. The 
terms corresponding to the harmonic oscillator operators of one-particle type may be 
easily obtained by means of the eigenvalue equation 


i (-V" + 4) Ui ) = (e + I) ) 


( 8 ) 
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Here 0e(C ) is the single-particle harmonic oscillator function in argument ( with os¬ 
cillator quantum number e. So, since many-particle oscillator shell model functions 
can be expressed in terms of single-particle harmonic oscillator functions, the following 
equality is rather straightforward 


^ ^ {{EKAjnry, ((e/j)A-i, {eij)A)rr'\\EKAjnT) 


(EKAjnT) 

J''T" 


{elj)A-i,{elj)A 

(ep)A-l’hb)A 


x{{EKAjnTyi{elj)y_,, 

X {{{elj)A-idelj)A)rT"\ 


{elj)y)rT"\\EK'A'jnT) 


— 

V(^) 





= E E 


(EKAjnT) hh)A-i,(eb)A 
J"T'' 


A^Ca-i + |) -|- A^Ca + |) 


X {{EKAJET); ((e/j)^_i, {elj)A)rT"\\EKAjnT) 


X {{EKAjnTy,i{elj)A-iAelj)A)rT"\\EK'A'jnT) . 


(9) 


The expectation value of the two-particle operators from Eq. (0) can be obtained 
by expressing them into the single-particle form. This may be accomplished by means 
of Jacobi coordinates 


^-1 = ■^{y(A-i) + ViA)) 

( 10 ) 

6 = -^{y(A-i) - y(A)) ■ 

With their use one can rewrite the operator as 


■(V(A-l) - V(A))^ + {.y{A-l) - y{A)Y — 2 I (“'^( 2 ) + 


( 11 ) 


The same transformation should be accomplished for the wave functions too. The 
antisymmetric two-particle oscillator shell model functions may be expanded in terms 
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of vector coupled products of the 'ipeij-Ktmjmt{^ 2 ) functions depending on the intrinsic 
Jacobi variable ^2 with the t/’(eO_i(^-i) functions depending on the Jacobi coordinate 
with the nonpositive index 0 


|((e/j)A-i,(e/j)A)J"T") = ^ |((e0-i,e/j7rf) J"T") 

(el) — i,elj 7 rt 

X (e/j)A)J"T"||((e/)_i,e/j7rf)J"T") . (12) 

Here eljntmjrrit characterize the separated two-particle subsystem, and Jacobi coor¬ 
dinate with nonpositive index is proportional to the centre-of-mass coordinate of this 
separated two-particle subsystem. Now, the functions 'il 2 eij-Ktmjmt{^ 2 ) are eigenfunctions 
of the operator 1 {-%> ^( 2 )) well cLS orthogoncil cinci normciliz6ci, thus for tho 

two-particle term of Eq. 0 one has 


l((e/j)A-n 

= E 


{elj)A)J"T"\ (H - l)i [-(V(A-l) - V(A))" + (^(A-l) - f(A))' 
(((e/j)A-i, (e/j)A) J"T"|((e/)_i, eljnt)J"T") 


{eiy_j^,{eljwty 


X mi)'.,-,, {elj)\)rr'\{(el)A, {elj7,tY)J"T") 

X 2(A - l){((ei)-,, elj^t)J"T"\^ (-Vfj, + (,%) |((ei)'-,. (eljntY)J"T") 

(el) — i,elj 7 rt 

X (((e/j)Ei, {eljyA)rr'\{{el).,,elj7it)rr') . (13) 


Here e is the number of oscillator quanta corresponding to the intrinsic motion of the 
separated two-particle subsystem. It is convenient to extract the coefficient 
2 (a 1 — 1) out of the Eq. ([I^) and move this coefficient to the right-hand side 

of Eq. (^. After that it obtains the form 
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X! X! [A{eA-i + Ba) + 3] 

[EKAjnT) (elj)A-i,{elj)A 
J"T" 

X {{EKAJUT); ((e/j)^-!, (e/j)^) J"T"|JilT) 


X {(EKAJET); ((e/j)^-!, (e/j)^) J"T"|JilT) . (14) 


Taking into account Eqs. 
element becomes 


([l^) and ([IlD the centre-of-mass Hamiltonian matrix 


{EKAjnT\^^§AE\EK'A'jnT) 

hui 


= E {{EKAjnT);i{dj)A-i,{elj)A)rr'\\EKAjnT) 


(EKAJHT) 

J”T" 


(elj)A-i,(>^lj)A 

(eb)h-l’(®b)h 


X {{EKAJET); ((e/j)A_i, (e/j)'^) J"T"| JilT) 


X {[2 i^A-i + Ca) + 2)]Eb)A-i,(eb)h-iEb)A,(eb)h 

-(H-1) ^ e (((e/j)A-i,(e/j)A)J"T"|((e/)_i,e/j7rt)J"T") 

elj 7 rt,{el)-i 

X (((rfj)'^_i,(eij)A^"T"|((d)-i.djA«)^"r")} ■ (15) 

It is possible to express the sum over the single-particle energy in terms of the 
total energy E. For this, let us express the total Hamiltonian, dehned as the sum of 
single-particle operators, in terms of the sum of two-particle operators 


E 

hoj 


El + yfi)} - E%) 


A-i E (^h) + ^(i)) • 

*J=1 

i<j 


( 16 ) 


Now, the total Hamiltonian matrix element between the oscillator shell model functions 
assumes the form 


{EKAJET\ 


E 

hui 


EK'A'JET) 
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= ^i{t){EKAjnT\{h^A-i) + h^A))\EK'A'jnT) 

= X! X! f (Ca-i + I + Ca + f) 

(EKAjnT) (eO')A-i,(eO')A 
J"T" 

X {{EKAJUT); ((e/j)^-!, (e/j)^) J"T"|JilT) 

X {{EKAjnry, {{elj)A-i, {elj)A)J"T"\\EK'A^jnT) . (17) 

The oscillator shell model functions are the eigenfunctions of the total Hamiltonian, 
so the corresponding matrix element is by definition 


{EKAjnT\ —lEK'A'JHT) = (e + \A) 5 ka,k'A' ■ (18) 

Tiuj ^ ' 

Comparing Eqs. (pTl) and (IT^) one obtains 


E- 6 kA,K'A'— X! X! ^(Ca-i + Ca) 

(EKAjnT) Cb)A-i,(eZi)A 
J"T" 

X {{EKAJET); ((e/j)^-!, {elj)A)J"T"\\EKAjnT) 


X {{EKAJETy {{elj)A-i, {el])A)J''T"\\EK'A'JET) . (19) 


Inserting in the equality (^) the expression (p!9D and introducing the coefficients 
of expansion of the oscillator shell model functions in terms of the ones with singled 
out dependence on the intrinsic coordinates of two last particles (SCFP’s), we get the 
final expression for the centre-of-mass Hamiltonian matrix element 


{EKAJET\^^AA.\ek'^JET) = iE + f) Skak'A' - (H - 1) V e 
X ^ ^ {{EKAJETy ((e/)_i, el]'Kt)J''T''\ \EKAJET) 

(eZ)_i, J"T" JeKAIWe) 

X {{EKAJETy ((e/)_i, elj'nt)J"T"\\EK' A'JET) . 


( 20 ) 






















3 Eigenvector problem of the real symmetric ma¬ 


trix 

For the centre-of-mass Hamiltonian matrix obtained, the general eigenvalue problem 
should be solved. Usually, this is carried out by some numerical methods. However, nu¬ 
merical methods are inefficient for highly degenerate subspaces, what is the case. The 
problem at hand is really less complicated, since the eigenvalues of the matrix under 
consideration are known. Moreover, this problem can be generalized, since eigenvalues 
of the quadratic Casimir operators for almost all physical interesting groups are well 
dehned. So, the problem can be formulated as to obtain the eigenvectors of the real 
symmetric matrix M, when its eigenvalues are known. We present here a new method 
which allows one to hnd the analytical expressions of the orthonormalyzed system of 
eigenvectors, corresponding to the chosen eigenvalue of the matrix. Since the sym¬ 
metrical properties of functions are usually formulated in terms of eigenvalues of the 
quadratic Casimir operators, this method makes it possible to obtain the set of func¬ 
tions with necessary symmetrical characteristics. Let us outline shortly the essence of 
the method. 

It is well-known that the real symmetric matrix M is expressible as the spectral 
decomposition of idempotent matrices [§] 

M = ^ . (21) 

a=l 

Here (Pa)nxn is the idempotent matrix of the rank Tq,, and Aq, denote the distinct 
eigenvalues of the matrix M. Idempotent matrices are known to obey the condition 

Pa ■ P/3 = f^a/sPa • (22) 

Using the Cayley-Hamilton theorem |^, the matrix P^ may be expressed in terms 
of the original matrix M and its eigenvalues 
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( 23 ) 


p» = n(v-M)/n(^/s-^«)' 

«=i «=i 

It is known, that the spectral decomposition of an idempotent matrix is [§] 

P„ = F„ • F+ , (24) 

where the matrix F^ is composed of eigenvectors of the matrix corresponding to its 
eigenvalues equal to one. In the case of real symmetric matrices, F+ is the transposed 
matrix. The eigenvectors satisfy the usual normalization condition F+ ■ F„ = 1. 

Taking into account the above spectral decomposition of the matrix it is straight¬ 

forward to show that 


M ■ F„ = F„ ■ . (25) 

It means that columns of the matrix Fq, are eigenvectors of the matrix M corresponding 
to the eigenvalue Aq,. 

The problem of constructing the eigenvectors for an idempotent matrix was once 
resolved in the case of coefficients of fractional parentage for antisymmetric wave func¬ 
tions 0. The method is based on the observation that the spectral decomposition of 
the matrix P is defined only within an orthogonal transformation. So, the matrix F 
defined by 


F = F • G , (26) 

also gives the spectral decomposition of the matrix P. The orthogonal matrix G has 
r(r — l)/2 independent parameters, they may be chosen in a way which allows one to 
fix the corresponding number of the matrix F elements. The best choice is 

{Fo)ij = 0 if 1 < i . (27) 
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In such a case the expressions for eigenvectors Fq, of the real symmetric matrix M 
corresponding to the eigenvalue Aq, take the form 


= 

(-Pa)ll ) 


{Fa)jl = 

iPa)lj 

iF^)u 


Wi = 

{Pa)ii - 




k=l 

{Fa)ji = 

{Fa,)a \ 

{Pa)ij '^{Fa)ik{Fa)jk\ 
k=l J 

1 = 2,3,.. 

., r and the corresponding set of j 


(28) 


Here n is the dimension and r is the rank of the matrix Pq,. The r columns of the 
matrix Fq may be obtained from the set of any r linearly independent rows of the 
matrix Pq. Consider the hrst nonzero row of the matrix Pq. Let it be the fc-th. 
Then the calculation begins with this fc-th row of the matrix P^: (Fa)n = {Pa)kk and 
{Pa)ji = Because the matrix Pq, is symmetric, every next i-th row calculation 

starts from the diagonal element (Fq)m. This Tth row of the matrix P^ will be linearly 
independent of rows just calculated if the obtained diagonal element {Fa)ii ^ 0. In this 
case the corresponding column of the matrix Fq, could be obtained. If this condition 
is not fulhlled, the next row of the matrix Pq should be tested. Since the matrix Pq 
has rank r, it will surely contain r linearly independent rows. Positive values of Fa are 
convenient, because the overall sign of the eigenvector is arbitrary. 

The outlined method for constructing the eigenvectors with necessary characteris¬ 
tics can be applied for diagonalization of the matrix. The obtained eigenvectors 

F will be the sets of coefficients a^'^oor introduced in Ref. 0. Then the coefficients 
of expansion of the oscillator shell model functions in terms of H-particle oscillator 
functions with singled out dependence on intrinsic coordinates of two last particles and 
with eliminated spurious states (CESO’s), as well as the intrinsic density matrices can 
be calculated. 
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4 Numerical benchmarks 


Let us use as a benchmark the formation of the matrix ( ), its diagonalization, 

and calculation of the CESO’s and intrinsic density matrices for the A = ?> case. 

The eigenvalues of the matrix makes the sequence |, | • • • E —EminT |, where 

E is the total oscillator energy value of the state. 

Let us start from the case E = Emin. Then, for three particles there is only one 
conhguration Ki, and the total angular momentum as well as the isospin can take only 
one value 

A-. ^ (ooi)’: JUT = Vi , 

The state (1) in the conhguration space will be denoted by the sequence of sing¬ 
le shells (e/j)"' accompanied with JT values in correspondence with the momentum 
coupling scheme. 

(1) ((ooq“)iy 

In the case when E = Emin, the J and T could take only the displayed values. 
For brevity, let us deal below with the same JUT = values also for higher total 
oscillator energies. 

In the case of three particles and E = Emin = 0 the matrix of the operator Ef^ ^ 
is the single number 


/ -^c.m. \ _ 3 

\ hu / A=3,£;=0,Ji7T=i + i ^ 

This is why the oscillator shell model states characterized by Emin contain the centre 
of mass in its ground state and all r — 1- 

Now the centre of mass is not excited, hence the CESO’s coincide with the 
coefficients of expansion of the oscillator shell model function in terms of the ones 
with singled out dependence on intrinsic coordinates of two last particles. They are 
presented in Table 1. Here the more usual spectroscopic notation is used to 

denote the state of a subsystem containing two separated particles. 

The intrinsic density matrices, in the case when H = 3, E = 0, JT = are 
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W('s„)=(i) , W(='Si) = (i) . 

In the case oi E = -Emin + 1 = 1, there are two conhgurations and a number of total 
angular momentum and isospin values appear: 


A-. ^ (ooi)hiih^ 
A'2 ^ (ool)hlli): 


JET 


1-1 1-3 3-1 

2 2’ 2 2’ 2 2’ 


JET 


1-1 3-1 3-3 5-1 
2 2’ 2 2’ 2 2’ 2 2‘ 


Now, the JET = I"*"! values may be obtained for the three states in the conhgu- 
ration space: 

( 1 ) ((ooybl,(lli)ii)ii, 

(2) ((00i)h0,(lli)ii)ii, 

(3) ((00i)h0,(lli)fi)ii. 

In order to calculate the ( E^ ^ ) matrix, the coefficients of expansion of the oscil¬ 
lator shell model functions in terms of the ones with singled out dependence on intrinsic 
coordinates of two last particles for E = -Emin +1 should be known. They are presented 
in Table 2. Here only the values of the indices different from the preceding row are 
shown. 

Since there are three states in the conhguration space, the dimension of the matrix 
of the operator E^ ^ will be three: 


E. 


c.m. 


huj 


A=3,E=l,jnT=^ i 


I 2 


V 


3J2 \ 


1 _2_ 

6 3+2 

U 

9 9+2 

18 '' 


where the rows and columns of the matrix ( E^. ) are labeled by the states in the 

conhguration space that are taken in the sequence (1), (2), (3). Due to the symmetric 
nature of this matrix here we show only its upper triangle. The rank of the idempotent 
matrices of this matrix is two, so there will be only two normalized eigenvectors. After 
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diagonalization of the matrix ( ) according to the outlined above method and 

using Eq. (^8[) , one may obtain 


[ 7^2 0 \ 



1 4 

3V2 3V2 



where rows of the matrix ( aKAflo,r ) are labeled by the states in the configuration 
space that are taken in the same sequence (1), (2), (3) as for ( ), columns are 

labeled by the additional integer quantum number F which can take two values: 1 and 
2, in correspondence with the rank of the matrix ( ) or, in another words, with 

the number of its linearly independent eigenvectors. It should be stressed, that the 
empty place instead of the hrst element of the second column means that its value is 
equal to zero, in full accordance with prescriptions of the method. 

The results of CESO’s calculations are presented in Table 3. It should be stressed, 
that CESO’s are labeled with the set of indices: {EKAJF[T){{el)_i,elj7rt)J"T'' and 
the additional integer quantum number F. There are only different sets of them, but 
they can belong to all states in the conhguration space, regardless of the origin of the 
ground conhguration. For example, if seven different sets of CESO’s indices are taken 
from the state (1), then only one can be taken from the state (2), since the remaining 
six have the same sets of indices as in the state (1). At last, a set of indices will be 
taken from the state (3), since the remaining four may be met in states (1) and (3). 
Table 3 shows only the values of the indices different from the preceding row, whereas 
CESO’s which are absent here have zero values. 

Finally, by means of obtained CESO’s the intrinsic density matrices can be calcu¬ 
lated. 

The intrinsic density matrices in the case of A = 3, E = 1, JT = 
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Wf"Pn) = 


12 


V 


, W('Pi) = 


- I 1 

4 


, W(^Pi) = 


V 


The rows and columns of the intrinsic density matrices are labeled by the additional 
integer quantum number P, preserving ascending order, so the dimensions of the ma¬ 
trices can vary from 1 to 2. Only the upper triangle of the intrinsic density matrices is 
shown due to its symmetry under index permutation. 

In the case of P = Pmin -l- 2 = 2 there are six conhgurations that produce a number 
of total angular momentum and isospin values: 


A'. ^ (00i)h20i): 

JUT 

A, ^ (00l)h22f): 

JHT 

A, ^ (00i)h22|): 

JHT 

A. ^ (ool)'(lli)t 

jnr 

A, ^ (ooi) (lli) (llf): 

JHT 

/Q ^ (ool)' (lli)': 

JHT 


1+1 1+3 3+1 

2 2 ’ 2 2 ’ 2 2 ’ 

1+i 3+1 3+3 5+1 

2 2 ’ 2 2 ’ 2 2 ’ 2 2 ’ 

3+1 5+1 5+3 7+1 

2 2 ’ 2 2 ’ 2 2 ’ 2 2 ’ 

1+1 1+3 3+1 

2 2 ’ 2 2 ’ 2 2 ’ 

1+ 1 1 + 3 3+1 3 + 3 5+ 1 5 + 3 

2 2 ’ 2 2 ’ 2 2 ’ 2 2 ’ 2 2 ’ 2 2 ’ 

1+ 1 1 + 3 3+1 3 + 3 5+ 1 5 + 3 7+1 

2 2 ’ 2 2 ’ 2 2 ’ 2 2 ’ 2 2 ’ 2 2 ’ 2 2 


Here and bellow the remarks, about notations mentioned above, will be valid, so 
the data will be presented in more concise manner. 


The JITT = i+i 


values could be obtained for nine states in the configuration space: 


(1) ((00i)bl,(20i)ii)ii, 

(2) ((00i)h0,(20i)ii)ii 

(3) ((ooyho,(22|)|i)ii, 

(4) (K)M.(ib)bi)M. 

(5) ((00l)ll,(lli)h0)ii, 

(6) (((ooi)M.(iihM)io.(iii)lhb. 

(7) (((ooi)M.(ii|)M)ii.(iii)lhH. 

(8) ((00i)ii,(llf)bl)ii, 

(9) (K)M.(iii)bo)M- 


15 









The number of the coefficients of expansion of the oscillator shell model functions 
in terms of the ones with singled out dependence on the intrinsic coordinates of two 
last particles is 92, so they will not be presented here for the reasons of space. 

Since there are nine states in the conhguration space, the dimension of the matrix 
of the operator will be nine 


hui 


jnT=\+\ 

A=3,E=2 


1 

1 

2 

2 

1 

5 

'svT 

3\/6 

3vT 

3^2 

3^3 

3v45 

1 

1 

2 

2 

1 

5 


9vT 

9vT 

9\/2 


9v45 

0 

20 

5 

5 

0 

_ 2 

9v^ 

9v4S0 

9vT0 

9v4s 

7 

3 

1 

6 

0 

_ 2 
3^3 

0 

0 
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9 

2 

9 

0 

0 

0 



7 

_1_ 

1 

5 



3 

6^3 

3^2 

9vT0 




2 

_ 1 

5 





3v4; 

3v^ 





7 

_ 5 





3 

6vT 
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The rank of the idempotent matrices of this matrix is 4. The diagonalization of the 
matrix ( ^ ) is accomplished within the framework of the outlined above method 

using 


jnT= 


1 + 1 


j \.o 11 ± 2 2 

\ O‘KA;00,r )a=3,E=2 — 


^ V2 

_]_ 

3\/2 

0 

0 
2 

9\/3 
4 

9\/3 

_4 

9 


V- 


0 

20 

9v^ 


2 

3 

0 

_ ]_ 

vT 

1 

9vT 

2 

9vT 

_ 2 _ 

9^2 

_ ]_ 

5 

9vT5 


0 

1 

0 


9vT0 

5 


9vT0 

5 

3^1) 

0 

2 

9 


0 

0 


1 

3 

_ 1 
3 

1 

3 

1 

1 

3\/2 

5 

’3vT0 
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Here the primes differentiate the notations of the states with coinciding orbital mo¬ 
menta. 


So far, for all excitation energies under consideration, only the subspaces corre¬ 
sponding to the nonspurious motion were extracted. It should be mentioned, that the 
same procedure can be carried out also for other subspaces corresponding to the pos¬ 
sible excitements of the centre of mass. For this, only the value of excitation energy of 
the centre of mass (the one from the sequence |, |... i? — i?min + |) should be changed. 
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5 Conclusions 


The proposed new method for elimination of ’spurious’ states enables one to produce 
two-particle oscillator intrinsic density matrices without exploiting the ICFP’s which 
untill now are not available for an arbitrary number of oscillator quanta. The distinct 
advantage of this method is in the complete rejection of group-theoretical classihcation 
of antisymmetric and translationally invariant many-particle states. So, a beneht could 
be gained due to simplicity and comprehensibility of such kind calculations. Another 
distinct feature of the method is that it do not involve any numerical diagonalization 
and orthogonalization, i.e. it is stable numericaly. So, the scale of calculations is 
determined only by available computation capabilities. Software tools based on a root 
rational fraction expression of numbers for calculation of CESO’s and two-particle 
oscillator intrinsic density matrices were developed and implemented in a computer 
code. The theoretical formulation have been illustrated by calculation of two-particle 
oscillator intrinsic density matrices in the case of A = 3, E = Emin, Emin + 1, E min -l- 2, 
and jnr = 

The presented procedure for projecting off the required oscillator shell model sub¬ 
space can be generalized by proposing a new method for calculation of analytical ex¬ 
pressions of eigenvectors if the analytical expressions of eigenvalues and entries of the 
real symmetric matrix under consideration are known. Since it is the case for all 
quadratic Casimir operators, this method can be applied to produce large irreducible 
subspaces with desired transformational properties. 

The new method for elimination of ’spurious’ states proposed in this paper enables 
one to formulate the problem of calculation of two-particle oscillator intrinsic density 
matrices in such form where this problem is no more complicated as matrix algebra 
and hence can be solved in a simple and efficient way. 
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Table 1. 


The coefficients of expansion of the oscillator shell model functions in terms of the 
A-particle oscillator functions with singled out dependence on the intrinsic coordinates 
of two last particles and with eliminated spurious states for the three-nucleon system: 

B = 0, A-=(00i)hi A = l, jnT={*\. 
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Table 2. 

The coefficients of expansion of the oscillator shell model functions in terms of the 
ones with singled out dependence on the intrinsic coordinates of two last particles for 
the three-nucleon system; E = 1, K = 01 ^ = 1, JHT = | ^ • 
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Table 2 (continued). 


The coefficients of expansion of the oscillator shell model functions in terms of the 
ones with singled out dependence on the intrinsic coordinates of two last particles for 
the three-nucleon system: E = 1, K = 10 ^ = 1, JUT = \ \- 
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Table 2 (continued). 

The coefficients of expansion of the oscillator shell model functions in terms of the 
ones with singled out dependence on the intrinsic coordinates of two last particles for 
the three-nucleon system: E = 1, K = 10 (nf) ||, ^ = 1, JUT = | |- 
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Table 3. 


The coefficients of expansion of the oscillator shell model functions in terms of A- 
particle oscillator functions with singled out dependence on the intrinsic coordinates 
of two last particles and with eliminated spurious states for the three-nucleon system : 
A = 3,E = 1, JT = II . 
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Table 4. 


The coefficients of expansion of the oscillator shell model functions in trems of A- 
particle oscillator functions with singled out dependence on the intrinsic coordinates 
of two last particles and with eliminated spurious states for the three-nucleon system : 
A = 3,E = 2, JT = ii . 
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Table 4 (continued). 


The coefficients of expansion of the oscillator shell model functions in trems of A- 
particle oscillator functions with singled out dependence on the intrinsic coordinates 
of two last particles and with eliminated spurious states for the three-nucleon system : 
A = 3,E = 2, JT = ii . 
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